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I. INTRODUCTION 

Recently there has been a discussion in the literature about when an imaging experiment carried out with spon- 
taneous parametric down conversion (SPDC) can be simulated classically [1,2]. We wish to discuss this problem in 
the general setting of quantum information theory and to give a class of experiments in which the measurements 
with a given entangled source can be reproduced using a classically correlated source. In addition, we show that 
certain aspects of the "ghost" interference experiment [3] can be reproduced using a classically correlated source. It 
is very important to have a clear understanding of what is being simulated. We show that while certain results can 
be reproduced by classically correlated sources, those results that depend on the phase of the two photon amplitude 
or biphoton can not. For a review of quantum imaging see ref. [4,5]. 

There has been a great deal of work done on determining whether a given state is entangled. The problem discussed 
here is different, we ask when the results of the experiment using a source that produces an entangled state can be 
reproduced by replacing the entangled state by a separable state. In the imaging experiments an entangled initial 
state is generated by SPDC and is transformed into the output state by the fields propagating through a linear optical 
system. The most important point is that the system has loss in it, indeed, the image is formed by removing photons 
from one of the beams. Loss may transform the state from an entangled state to a classically correlated state. When 
this happens all the measurements can be reproduced by a classically correlated source. However, this is not the 
general case since, even after losses, the state may remain entangled. Consequently, only certain measurements may 
be reproduced classically. This is the case we wish to discuss in this paper. The problem of substituting a stochastic 
classical source for a quantum source in optics experiments has been discussed many times. In the context of this 
paper, two reviews by Belinskii and Klyshko [6] are of interest. 

It is important to be clear what we mean by an experiment. The original imaging experiment using SPDC [7] 
consisted of a source of entangled photon pairs with each photon passing through a linear optical system to detectors. 
A measurement consisted of the coincident detection of the pair. Each photon was detected locally and the coincidence 
triggering corresponded to classical communication. The outcome was summed over many such pairs. In the original 
experiment one detector scanned the image plane but this can be eliminated by placing a CCD detector in the image 
plane. In this case the experiment is series of measurements performed on a fixed setup that measures the intensity 
correlation function (/^(ri)/^^ )) where the / 's are the normally ordered intensities and j = 1, ■ ■ ■ , N indexes the 
detectors in the CCD array. Note that these are N commuting operators. In the experiment of Bennink, et. al. [1] 
the source is replaced by a laser, beam splitters, a chopper and a moving mirror. In this case each measurement 
corresponds to a point by point scanning of the object plane. Light detected behind the object triggers the image 
plane detectors to their on state and the detector corresponding to the image point responds. The two experiments 
are related by replacing the SPDC source by a classical source. The classical source is meant to simulate a stochastic 
source that has a large bandwidth in the transverse wave numbers. If this bandwidth is the same as the bandwidth 
for SPDC, there will be no image behind the object plane. 

These experiments should be contrasted with a test of a Bell inequality [8]. Such a test consists of a series of 
experiments for which the source is fixed but the measurement operators are different. This is important because 
the quantities measured at each site are non-commuting operators, so they can not be measured simultaneously. It 
has been recognized for a long time that any experiment which only measures a fixed set of local operators can 
be reproduced by a local hidden variable theory in which the probabilities for each outcome of the experiment are 
themselves taken as the local variables. Then any correlations that exist can be modelled by separable density matrices, 
a simple example is given in section III. The experimental results of [1] are consistent with this conclusion. Some 
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SPDC experiments can be understood in terms of a picture first proposed by Klyshko [9,10] in which one detector is 
viewed as a source of photons propagating backwards in time, reflecting off the SPDC crystal, and then propagating 
forward in time to the second detector. It is often not appreciated that this picture has a rigorous mathematical 
foundation for maximally entangled photon pairs. This picture is genralizcd in VI C so it can be used to explain the 
experiment similar to that in ref. [1]. However, as can be seen, results that depend on fixed phase relations between 
two-photon amplitudes can not be simulated. Furthermore, it should be remembered that the quantum or classical 
nature of a state is determined by the statistics of the experiment, so even if the average of a set of measurements 
can be reproduced classically, it does not follow that the fluctuations in correlation measurements can be reproduced 
by classically correlated systems. In the imaging experiment there are two important results such that, as we shall 
see below, one can be simulated but the other can not. For a different perspective on this problem see ref. [11]. 

In the following we will analyze the statements made above and give examples of when an experiment using an 
entangled input state can be simulated by a classically correlated input state. We first review some well-known results 
about measurements and entanglement, and apply them to the simple case of two qubits to illustrate some of the 
issues discussed above. Then we discuss the imaging experiments. 



II. BACKGROUND 



For all the experiments that we discuss assume that the source produces a pair of subsystems A and B in a 
state pab • The subsystems are sent to two separate sites, also denoted by A and B, where local measurements are 
made. The Hilbert space describing the system is H = Ha <8> Hb- The set of observables is defined by {O^'b = 

O^ 1 £S> O^b <= S} acting on H where the sets of Hermitian operators {O^} and {O^} act on Ha and Hb, 

respectively. A state gab is said to be separable or classically correlated if it can be written in the form 

E(u) _ (it) \ 
VuO A ® a B I 1 ) 

u 

where {p u } is a probability distribution, i.e. a set of positive numbers that sums to 1, and a a (&b) is a density matrix 
acting on Ha (Hb)- An experiment will be defined to be classical (more precisely, classically correlated) if there exists 
a separable state such that for all G S 

trO% ] pab = trO { 2i ] (tab (2) 
We may actually be content to simulate one observable that is a linear combination of direct product observables, 

Oab=J2^°ab- (3) 

ij 

All of these are correlation experiments for which the experimentalists agree which pair of observables they measure 
for a given realization of the experiment. For example if the polarization of a pair of photons is measured then 
i = H or V corresponding to the measurement of a horizontal or vertical polarization in A, and similarly for the 
measurement in B, j = H or V. The four possible projection operators arc not only pairwise commuting but are 
locally pairwise commuting, consequently, the experiment can be performed using one fixed set-up in each laboratory. 
On the other hand if, for example, i = H,V,R and L where R and L refer to right and left hand polarized light, 
then the observables in A do not all commute and can not be measured simultaneously 1 . There may be more than 
one classically correlated state that can simulate a given experiment. If this is the case, the set of such states will 
form a convex set [12]. A simple illustrative example for qubits is given in section III where these measurements are 
discussed in a geometrical context. 

We illustrate how an experiment on an entangled state can be simulated using a separable state by a simple example. 
Consider an entangled pure state 'S, let Oab = Oa ®Ob be an observable, and let {(f) j} be the set of eigenvectors of 
Oa- Let {xj} be an orthonormal basis of Hb and consider the arbitrary state 



1 Of course it is possible to use a beam splitter to send half the photons to the pair of detectors that measures H and V and 
the other half to detectors that measure R and L but this does not change the argument. 
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N 

* = C 3k^j®Xk- (4) 

j,k=l 

The measurement of Oab gives 

tr(0 AB \^)m = ^2^(x 3 \0 B \x 3 }\c 3 \ 2 (5) 

3 

where 0a4>i = K4>i- It is not difficult to see that we get the same result for the classically correlated state 

<r = '52\<i>j}( < t>j\®\uj}(uj\ = ^Cjklxk)- (6) 

j k 

Therefore, any simple local measurement, that is, a measurement of a direct product of local operators, can be 
simulated by a separable measurement. 

Now consider a different observable obtained by a unitary transformation applied to Oa, 0' ab = 0' A ® Bl such 
that 0' A = WOaU and Oa do not commute. For example if subsystem A is a spin- 1/2 particle and Oa — &z take 
O'a = <J X . It is not difficult to show that the new measurements are not produced by the same classically correlated 
state since (<f>k\0' A \<f)j) ^ ^j$kj in general. Finally if {O a ' b , i,j € S} where the {O^} are pairwise commuting 
operators, this result immediately generalizes since the {<j>j} can be taken as simultaneous eigenvectors of all the . 

It is worth noting that if trOAB = we can always find a separable state cr(so) such that trOABP = 
(1 / 's )trO abc(s ) . To see this let p be the given entangled state on an TV dimensional Hilbert space, I N the identity 
operator, and D be the completely random state, D a = (l/N)Ijy [13]. The convex combination a(s) = (1 — s)D + sp 
is a state of the system which is entangled for s = 1 and classically correlated for s = 0. Therefore, there exists an 
s e (0, 1) such that the state a(s) = (1 — s)D + sp is classically correlated for s < s [14]. The reason for mentioning 
this is that many experiments do not look at absolute values but at relative values for which the normalization of the 
density matrices is not important. 

In a typical experiment, the state of the system measured, pab, is related to an input state, Pqab, by a completely 
positive, trace preserving, linear transformation 

Pab = T(p ab)- (7) 

The Kraus representation of T [15,16] is 



Pab = T(poab) = J2 V kP0A B Vl (8) 
k=i 

where the set of operators {14, k = 1, • • • , n} acting on H satisfy 

n 



fc=i 



with Ih is the identity operator on H. In the experiments of interest in this paper, each of the operators Vk is a direct 
product of operators acting on Ha and Hb, 

Vk = A k ®B k , (10) 
In quantum optics experiments the Vk can be computed from classical optics. Using eq. (8) we have 

n 

trO A BPAB = E tr °AB (VkPOABV^ 
k=l 

n 

= Y, tr ( V k°ABVk)p AB (11) 



k=l 



the first line is referred to as the Schrodinger picture and the second line as the Heisenberg picture. Using our locality 
assumptions 
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VlO AB V k = A\0 A A kl ® Bl 2 B B k2 k = {k u k 2 ). 



(12) 



We can now rewrite eq. (2) as 

trO^PAB = £ft, (trO^a^) (trOfa^) . (13) 

u 

If both O^* and are positive operators then the right hand side of eq. (13) is the sum of positive terms and 
there can be no relative phase information about the subsystems extracted from this expression. Then since any 
Hcrmitian operator can be written as a difference of two positive operators, there is a stringent condition on the type 
of experiment that can be simulate by classically correlated systems. 

It will be useful to record a few well-known results for pure states. Let ^ be a pure state with Schmidt representation 
given by 

N 

* = 5^c fc |0 fc )<8> \ Xk ) (14) 
fe=i 

[16]. Recall that the set {|cfc| 2 ,/c = 1,---,N} is a probability distribution and the states {4>k} and {xk} are 
orthonormal sets in Ha and Hb, respectively. The state ^ is separable if c k = Ski, otherwise it is entangled. If 
Ck = 1/y/N, the state is said to be maximally entangled. For a maximally entangled pure state vl/( M ) a Kraus 
operator applied at one subsystem can be replaced by an operator applied at the second subsystem, 

(A, ® B r ) * (M) - (1a ® B r A T S ) *< M ) , (15) 

where A s <pj — J2k <Pk [A s ) k j and A^Xk = J2j Xj (A s )kj ■ I n experiments it is often useful to place objects corresponding 
to Kraus operators such as beam splitters or lenses in both laboratories A and B, however, for maximally entangled 
states this is entirely equivalent putting them all in one of the laboratories. This fact is exploited in remote state 
preparation [17-19]. 



III. ENTANGLEMENT WITNESSES 



For finite dimensional systems the question of whether a given state, represented by an N x N density matrix, is 
entangled can be answered by operators, called entanglement witnesses [20,21], in the space of Hermitian operators. 
Entanglement witnesses have a geometric interpretation as hyperplanes that separate a given entangled state from 
the separable states [20-22]. In this geometric picture, an optimal entanglement witness is one which is tangent to 
the surface of the separable states. Recently, the optimal entanglement witnesses have been characterized in terms of 
a set of coordinated local measurements [23,24]. An entanglement witness, W, is defined for an entangled state p as 
a Hermitian operator such that 

tr(pW) < 0, (16) 

and for any separable density matrix t 

tr(rW) > 0. (17) 

The entanglement witness is said to be optimal if there exists a separable density r such that equality holds, that is, 
W determines a plane tangent to the separable states. 

Consider the simplest case of two qubits [23,24]. We shall denote the Pauli matrices by cr^, fj, = 0, • • • , 3 where er 
is the 2x2 identity matrix. Consider the state 

p=|$ + )(<D + | (18) 

where $ + = ( 1 00) + |11)) is one of the Bell states, |0) corresponds to the spin-up state, and |1) the spin-down state 
in the basis in which the Pauli matrix cr 3 is diagonal. It has been shown that the optimal entanglement witness for 
this state is 

W = | ((7 <8> cr - 3r ) (19) 
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where To can be written in terms of local projection operators, the factor of 2/3 is a normalization. For j — 1,2,3 
define the projections 

Pf = \ K ± aj) (20) 

and 

Pf = P? ® pj>, (21) 

then 

ro = l [(P 3 ++ + Pf~) + + P") + + P2 + )] • (22a) 

We see that To can be written as the sum of three pairwise commuting operators that are, however, locally non- 
commuting observables corresponding to the measurement each spin along one of three orthogonal axes. For a fixed 
experimental set-up, one of the pairs of operators in the parentheses can be measured. For example, the measurement 
of the state in eq. (18) along the 1-axis can be reproduced using the separable state 

r = P++pP++ + P~pP~ = \ (P x + <8> + Pf <8> Pf ) . (23) 

Similarly, the other two measurements can be reproduced if the input state is a (different) separable state. However 
the measurement trWp < can not be reproduced by a separable state, r, since for all such state trWr > 0. 

We can exploit this geometric picture further to show that there may be many observables that will satisfy eq. 
(13). Let t be the nearest separable state to p in the sense that 

\\p-T\\ = y/tr[(p-T)*]>\\p-T \\ (24) 

for all separable states r. Then there is an optimal entanglement witness for p through t . Suppose t has rank 
greater than one. Then the hyperplane of the entanglement witness intersects the separable densities in, at least, a 
line containing r . For any separable state r on the line 

tr[(p-T )(r-T O )} = 0. (25) 

Using the fact that the separable density matrices form a compact convex set, we can write 

r-r = J2 a ^A ) ®^B ) (26) 

u 

(u) (u) 

where the a u are real but not necessarily positive and the a A <X> a B are extreme points of the line. Therefore, t — tq is 
an observable that corresponds to local measurements with classical communications, consequently, when the system 
is in the entangled state p its measurement of this variable can be simulated using the separable state r . 



IV. SPDC EXPERIMENTS 



The general theory of the type of experiment shown in fig. 1 has been given in [25]. The source discussed in that 
paper was maximally entangled in frequency and transverse wave number. It was pointed out that the source po 
could be envisioned as entangled in transverse modes which are not necessarily plane wave modes. In the present 
context these modes are generated by the POVM's acting on the maximally entangled transverse modes of the beams 
generated by the SPDC crystal. Specifically, suppose the pump is a plane wave of angular frequency lo p normally 
incident on a crystal that is infinite in the transverse direction. The state of the down converted biphoton is 

*o = fA^h P>lK - Wj, -Pj)|0) (27) 

3 

where a\{uj, p) is the creation operator for a photon at the output face of the crystal with polarization r, frequency 
w, and transverse wave vector p. Then we can write the state propagated to the detectors located at ri and r 2 



5 



*ab = ^2gA(ri,Uj,Pj)gB(r 2 ,u) p - Uj, -Pj)fjal(wj, Pj>!( w p - -Pj)l°> 
i 

= fj a l ^ 1 ' ^ ' P ^ ' - W 3 » ~Pj) |0> ( 28 ) 

3 

where the operators in second line are defined in the obvious way (cf. cq. (12)). For example, 

<4(ri,u;j,pj) = gA{r u Uj,pj)a\{wj,pj). (29) 

is the creation operator at detector one. The propagators g are determined by the linear optical systems in each 
arm of the experiment (see appendix VIA and [25] for detailed calculations). They can be expanded in terms of any 
convenient set of orthogonal transverse modes as is usually done when the field propagates in a fiber. It may be useful 
to recall that ^(ri, coj, Pj) is a transformation from the transverse momentum pj to the transverse position Y\t ■ 
The mode label j and the longitudinal variable t\l are parameters 

It is the entanglement of the pair in their transverse coordinates that is most important for image transfer. In 
many experiments the maximally entangled state is a good approximation. As shown in eq. (15), it is the maximal 
entanglement that makes it possible to put all the optics in one arm. If the pump is not a plane wave \fo is no longer 
maximally entangled in the transverse direction. 

In the quantum imaging papers using SPDC one produces an image or interference pattern in coincident counting. 
Since the SPDC has a large wave number bandwidth even for a narrow frequency range there is no image or interference 
pattern behind the object plane. Both of these effects can be simulated by classically correlated sources. A second 
important feature of the quantum case is the relation between the position of the aperture in one arm and the detector 
array in the second arm. These relations depend on the phase relation between the two-photon amplitudes or biphoton 
and, as we shall see, can not be reproduced by classically correlated sources. This is a specific example of the result 
discussed in connection with eq. (13). 

Starting with cq. (28) we can first ask what happens if we simply detect all the photons at B discarding all the 
photons at A 

PB = tr A \^AB){^AB\ 

= fjk a i( r 2,u p - Wj, -pj)|0)(0|ai(r 2 , (j p - u) k , -p fe ) (30) 

jk 



fjk = fjfk / d 2 r 1T g A {r 1 ,u v p J )g* A (r ll uj k ,p k ) (31) 



where 



and Tit is the transverse component of iv As pointed out in [1,2], if g A is unitary the fy. = \fj\ 2 5n- and the marginal 
density is the same as that of the source. We have not put indices on the g A and g B because we are concerned with a 
fixed experimental set-up. However, if we consider different set-ups, as we would in testing Bell inequalities, then we 
need to index them. 

We now consider coincidence detection. For simplicity we take our detectors to be polarization independent, point 
detectors oriented along the direction. The electric field operator for such a detector located at r is given by 

E = + £?(-) = £«t 

25(+>(r,t) ^Eja^p^"-^ k, = Pj +e L M - pj (32) 

3,1 

where Ej is determined by dimensional considerations, q is the polarization index, and pj << u/c. The probability 
of two detectors firing is proportional to 

tr (pabE^E^E^E^) = \A l2 \ 2 (33) 

where 

A 12 = (0\E[ +) Ei +) \^ AB ) = (0\EM( ri , h)EM (r 2 ,t 2 )|* AB ) 

= 9a(ti, Uj, Pj)gB(r 2 , oj p - Uj, -pAfjEjEhe-^'e-*"*** , (34) 
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u>h = Up — u)j, and = — pj + e^y — p 2 . An is the biphoton amplitude. 

In the imaging experiments one detector collects all the light that passes through the object, T in fig. 1, i.e. the 
correlation that is measured is the correlation when the transverse component of ri is integrated over. To carry out 
the computation we assume degenerate phase matching so that we can write lu± = uj p /2 + v and 0J2 = u p /2 + v' where 
frequency phase matching gives v 1 = —v. We require \v\ << u) p /2. In the paraxial approximation we can write 



A 



12 



J dvuivy^-^e^^+^l 2 ^ J d 2 xt(x)e- lK2 ^-^ 2 (35) 



where i(x) is the transmission function of T, and 



2c\ 2 / 1 1 1 



The term containing K comes from converting the sum over j into an integral in the amplitude A12, and is therefore 
an effect of the entanglement of the photons. When K — > 0, the Gaussian becomes a delta function and the coinci- 
dent counting rate becomes proportional to |f(x2)| 2 . The fact that the image and object planes are related by the 
geometrical optics condition K — was derived and experimentally studied in [10] where it was also generalized to 
the case of non-degenerate SPDC. 

We now ask whether this experiment can be simulated using a classically correlated source of coherent states. 
Consider the classically correlated state 

TAB =YlH T Al (37) 

3 

where {qj, j = 1, 2, • • •} is a probability distribution, 



AB - T A ii9T B T 



■%' = m')(W\ X = A,B (38) 



and 

m 
m 

The state \a m ) is a coherent state generated by the creation operator a\{w m , p m ) defined in eq. (27) and the coherent 
states \(3 m ) are generated by operators a\(oj p — u> m , — p m ). The restriction to coherent states is not strictly necessary 
to show that the results can be obtained from a classically correlated state, that is, we do not have to assume that 
the states sent to A and B are themselves " classical" . 

For the experiment illustrated in fig. 2, it is not difficult to show that 

C = T lim 1 £ dtr £ dt 2 tr (r AB E^E^E^E^) = ^ q 3 tr A (r^E^E^) tr B (r^E^E^) 

= Jd 2 pJ dV?(p,p')l.9A(ri,^ p /2,p)| 2 | ffB (r 2 ,^/2,p')| 2 (40) 

[1,2], where we have converted the discrete sum over modes to integrals and defined 

ErflV^| 2 |a m | 2 |/3 m f - dV¥<?(p,p'). (41) 
3 

Equation (40) is a special case of eq. (13) and contains no relative phase information between A and B. In arm B, 
the detector is placed in the focal plane of the second lens then, as shown in the appendix eq. (61), in the diffraction 
limited case 

| ffB (r 2 , Wp /2,p')| 2 =G^)( P '-glx 2 ) (42) 
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where x 2 is the transverse coordinate of one of the point detectors located at D2. Taking the correlation function 

g (p,p') = <5 (2) (£P'-p) (43) 



we get for the average number of coincident counts per pulse 

C = C |t(e^x 2 )| 2 . (44) 
h 

Note the detector in arm B and the object T are both in the focal plane of the lens. This can be interpreted as a 
point by point mapping of the mask onto the plane of the detectors at B. Equation (40) can be viewed as incoherent 
imaging [1,26]. The fact that a point by point imaging is at the heart of this type of experiment was noted some time 
ago in rcf. [10]. 

The classically correlated source modelled here is a stochastic source that is broadband in spatial frequencies. This 
condition is included to ensure that there is no image in the Dl plane, as can be seen by setting gs = 1 in cq. 
(40). The main characteristic of this source is that whenever it emits a beam to A with transverse wave number p 
it emits a beam to B with wave number p/e. The SPDC source has the same type of correlation but is significantly 
different because it imprints relative phase information on the beams in A and B which influences the correlation 
measurements. The important phase information exists in the transverse momentum dependent factors ip(p, X) in 
cq. (49) and eq. (50). These phase factors disappear in eq. (40). This places stricter conditions on where the images 
are located than in SPDC. 

The "ghost" interference experiment in ref. [3] displayed remarkable features. A simplified diagram is shown if fig. 
3. The first important feature is that an interference pattern can be seen in the coincidence counting but if single 
photons are detected in the plane of Dl no interference is seen because the SPDC is a broad band spatial frequency 
source. The second feature is that the period of the interference pattern is determined by the distance z from T back 
to the surface of the crystal and then to D2. For the two slit experiment the coincident counting rate was proportional 
to cos 2 (pd/2) where d is the slit separation and p = 2tt(x2/zX), x 2 being the position of a detector in array D2. A 
similar result holds for the single slit case. As the authors of ref. [3] point out: "The interference pattern is the same 
as one would observe on a screen in the lane of D2, if Dl is replaced by a point-like light source and the SPDC crystal 
by a reflecting mirror." This is just the Klyshko picture [9]. 

Let us examine this in a bit more detail. The detector Dl is at the origin of the plane at "infinity" so it detects 
the p = transverse mode. Since the SPDC is a broad band spatial source there are many different modes impinging 
on the slits. When a photon is detected by one of the detectors in D2 one of these modes p' is selected and we 
get coincidences when p' w 2irn/d where d is the distance between the slits and n is an integer. For a given p' the 
maximum of the interference pattern in the plane at Dl corresponds to w p' / (u>/c). Because of the broad bandwidth 
in p' , this interference pattern is averaged over and is not visible without the coincidence measurement. 

The interference pattern observed in [3] can also be simulated using a similar classically correlated source as is 
shown in VI B. As for the case discussed in [1], not all the features can be reproduced and the placement of the 
lenses LI and L2 in fig. 2 play a critical role. It is important to understand what about these experiments is being 
simulated. In both cases the pattern, image or interference fringes, and the fact that there is no image or interference 
if one looks behind the object are simulated. The results relating the placement of detector D2 and the object T can 
not be simulated because, as we saw above, because this arises from phase information stored in the entangled system 
and not separately in the subsystems A and B. As in the imaging case this phase information is that in the factors 
tp(p,X) in eq. (51) and eq. (52). From eq. (62) it is easy to see that for the double slit experiment the fringes are 
determined by p — 27r(x 2 //2A), that is the distance z in the SPDC experiment has been replaced by the focal length 
of the lens in arm B. 

In a recent paper [27] it has been shown that for an entangled source one can obtain both the image and its Fourier 
transform by only varying the position of a lens and the detector in the arm of the experiment that does not contain 
the aperture plane. This paper exploits the fact that for maximally entangled states one can locate the optics in 
either arm of the experiment. The authors explain their results by using entanglement, path indistinguishability, and 
complementarity in the transverse momentum and position. The complementarity argument must be made with care 
since it means different things to different people. The spatial and momentum (wave number) complementarity is 
a wave property and classical electromagnetic fields are waves. After all, one can display the spatial position and 
spatial frequencies of a source by simply changing the position of a lens. The path indistinguishability along with the 
quantum rule about adding amplitudes is more convincing since that invokes the particle properties of the photon 
which distinguish the quantum and classical pictures of the electromagnetic field. By using special cases of mixed 
state, the authors obtain a special case of eq. (13) and reach the same conclusions as those in this paper for the 
imaging case. However, they conclude, for their mixed states that they can not reproduce the interference fringes. 
This is not a general result, as we have shown. 
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We have shown that it is often possible to simulate some features of an experiment using an entangled source by 
a classically correlated source. Therefore, it is important when invoking entanglement as critical to an experimental 
result that the possibility of such a simulation is eliminated. Violation of Bell inequalities is a special case of using 
entanglement witnesses as observables to show that entanglement is necessary to explain some experimental mea- 
surements. Some quantum imaging experiments with SPDC sources are examples in which it is possible to simulate 
some of the measurements with classically correlated sources. In particular, the imaging pattern can be reproduced in 
coincidence counting even when there is no image or interference pattern visible behind the object. On the other hand, 
it is impossible to reproduce the aspects of the quantum experiments that depend on the relative phase information 
in the state of the entangled subsystems. This is a consequence of eq. (13). 

We have not answered the general question of what local observables can be simulated by a classically correlated 
source for a given entangled source. Finally, it should be noted that the converse problem of when a classically 
correlated state can be used to study entanglement has been studied extensively starting with ref. [29]. 
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VI. APPENDIX 



A. Notation 



We briefly review the general form of a coincidence counting experiment [25]. The positive frequency part of the 
electric field at time t at the input of a detector located at r =ze z + x is given by 



4 +) (r,t) = J J dutfpe-^EMgfru^afi (w,p), 



where ap (ui, p) is the annihilation operator at the source for a photon of angular frequency w, transverse wave number 
p, and polarization (3. The unit vector e z is the inward normal to the surface of the detector, and x is the transverse 
coordinate on the face of the detector. E{oj) is required for dimensional reasons and is a slowly varying function 
that can be taken to be constant in the experiments discussed in this paper. Finally, g (r, u, p) is the optical transfer 
function or Green's function determined by the classical linear optical system between the source and the detector. 
It may be worth noting that classically this quantity connects electromagnetic fields that live in real space-time while 
in the quantum mechanical context it is a mapping of operators that operate in Fock space. This dual nature is 
the source of many results in this paper and is also the source of confusion about the " photon" . It is important to 
recognize that linear superposition in classical electromagnetic theory plays an essential role in determining g. In the 
quantum mechanical case g is best thought of as being determined by the boundary conditions defining the modes of 
the system independently of the state of the system. That is, the modes are the same whether the system is in a one 
photon state or in a "classical" state for which modes are excited by coherent states. 

As an example of the computation of g, we consider arm A of fig. 1. We use the notation i(x) for the transfer 
function of the object where x is the transverse spatial vector in the plane of T. The thin lens transformation is 

J(x,/)=V(|x|,-~) (45) 

where / is the focal length and x is the transverse coordinate in the plane of the lens. Then, after integrating over 
the plane of Dl, 

&4(ri,u;,p) =Jd 2 x a Jd 2 x f Jd 2 x s t(x a )h UJ (x a -x f7 d' 1 )x (46) 
i(x / ,/)fc w (x / -x a ,di)e ip - x '. 
where r\=(d\ + d'{)e z , x s is in the output plane of the source, and in the Fresnel approximation 
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f —iuj\ e*( w / c ) d h) 1 

h ^ d) ={2^)— *<N£i> < 47 » 

iP(x,^P) = e l ? ( ^ Px \ (48) 
This notation is very convenient [25,28]. Evaluating g& gives 

ffA(ri,w,p) - (-*27r^) ^-^- e ^(*+^)^(p,-£D) x 

y d 2 x a t( Xa )e^f^t'- d '^^ a ,-kj^j-), (49) 

where D = di + . 

For arm B taking r2 = x + d 2 e z and performing the integration over x s assuming that source has a very large 
cross-section, we have for each plane wave input mode 

3B (r 2 ,^p) - J d 2 x s ^(x-x s ,d 2 )e^ = e^/ c ) d V p ->(|p|,-^ 2 ) (50) 

which is the Fresnel approximation to the plane wave exp(iy/ (^) 2 — p 2 <i 2 + p ■ x). Now consider what happens if we 
introduce an aperture function in this expression. 

B. Ghost interference 

Figure 3 shows a schematic of the ghost interference experiment [3] . In this case 

9a ((di + d'-^e^u)) = J d 2 x a j d 2 x s h 0J (-x a ,d' 1 )t(x a )h 0J (x a -x s ,di)e lp ' :XL3 . 

:J j7 J i^M^l,-f i )t(^)e-'"-t(\p[,--d,) (51) 



-( — 

\ 2-kc 



and 

g B (d 2 e z , J, p') = e^ w '/c)d2 e V-x 2V , ( | p /| 5 _ £ ^ (52) 



In the degenerate case, lo — u' may be replaced by uj p /2 in the ^-function and p' = —p. Then defining the biphoton 
wave function to be 

* = j 'd 2 p j 'duu(u)at{^+u,p)a\(^--u,-p)\0), (53) 



the amplitude is given by 

g-iw p (Ti+T2) 12 



A 
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J dvu{v)e MTl - T2) A (54) 



d[ 

A = A "f ^«V'(|x a |,g^-Mx M|x a -X 2 |,g- r J- r ) 



di + d 2 ' 



= Ao f d 2 x a ^(x a , ^ ( — L- + 1) )t(x Q ) e - 4 ^/ 2c )( X2 / dl + d2 )- x °V(|x 2 |, ( 55 ) 
J 2c \di+d 2 dj/ 2c c(i + d 2 

In the far field Fraunhofer approximation we may replace the ip's in the last line by 1 

A " = ^ AT*' < 56 » 

where t(p) is the Fourier transform of the aperture transfer function. 

In the classical source case we take arm A to be the same as in fig. 3. Then from eq. (51) we have 
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\g A ((d 1 +d[)e z ,u;,p)\ 2 =G A \ J d 2 x a ^(\ Xa \,^^-)t( Xa )e^\ 2 (57) 

In arm B of the experiment we use the set up shown in fig. 2, 

5s( r 2,w,p) = y y d 2 x ; d 2 x ;s ^(x2-x z ,/2)P(x z )/(x ( ,/ 2 )/i w (x z -x s ,rf2)e 4p ' Xs (58) 

where we have include a finite aperture for the lens, P{x{). To make the evaluation of the integral simple, we shall 
take 

P(x,) = e-l-'l 2 / 2 - 4 . (59) 

Then 

Mr 2 ,c,p)| 2 = ^1) A 2 e"( p -^-) 2A . (60) 

In the limit of large A, the exponential will be very small unless (p— ^-^-x 2 ) 2 « 1/A. In this case we may write to 
a good approximation 

\g B (r 2 ,u,p')\ 2 = f^l) A7r5 2 (p'-^|x 2 ). (61) 

Combining this result with eq. (57) in the far field limit, d\ — > oo in the integral, and the classical correlation function 
cq. (43) gives for cq. (40) 

C = G#(e^)| 2 (62) 

where the lens LI in fig. 2 is not necessary. 

C. Klyshko picture 

Consider an experiment like that shown in fig. 1. For a biphoton source the amplitude for the process can be 
written as 

A = j d 2 pg A (n , w, p)g B (r 2 , w, -p). (63) 

The Klyshko picture is based on the observation that .gs(r 2 , w, — p) =.9s(r 2 , — w, p)*, see for example eq. (52). The 
complex conjugate Green's function can be viewed as a source located at r 2 progagating backwards in time reaching 
the SPDC crystal with transverse momentum p. At the crystal the beam propagates forward in time to ri. 

To generalize this picture, suppose that the beams with transverse momentum p have independent arbitrary phases 
so eq. (63) now takes the form 

A = Jd 2 pg A (r u u,p)e ieA ^g B (r 2 ,uj,-p)e ie »(-*>\ (64) 

where 

( 6 *Mp)) = o ^ec(p) e -t0c( P '^ = KgWfp _ p ') C = A or B (65) 

( e *0A(p) e -»0B(p')) = 

and K is a constant. With these assumptions 

(\A\ 2 )=K J d 2 p\g A {r u u,p)\ 2 \g B {r 2 ,uj,-p)\ 2 (66) 
which is of the form of eq. (40) for the classical correlation q(p, p') = K5^ 2 \p' + p). 



11 



[1] R. S. Bennink, S. J. Bentley, and R. W. Boyd, Phys. Rev. Lett. 89, 113601-1 (2002). 

[2] A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko, and M. C. Teich, Phys. Rev. Lett. 87, 123602-1 (2001). 
[3] D. V. Strekalov, A. V. Sergienko, D. N. Klyshko, Y. H. Shih, Phys. Rev. Lett. 74, 3600 (1995). 
[4] Y. H. Shih, J. of Mod. Opt. 49,2275 (2002), 

[5] L. A. Lugiato, A. Gatti, and E. Brambilla, J. Opt. B: Quantum Semiclass. Opt. 4, S176 (2002) and quant-ph/0203046. 
[6] A. V. Belinskii and D. N. Klyshko, Laser Phys. 2, 112 (1992), A. V. Belinskii and D. N. Klyshko, Phys. Uspenkhi 36, 653 
(1993). 

[7] T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. Sergienko, Phys. Rev. A 52, R3429 (1995). 
[8] J. Bell, Speakable and Unspeakable in Quantum Mechanics , (Cambridge Univ. Press, London and New York, 1993). 
[9] D. N. Klyshko, Photon and Non-linear Optics (Gordon and Breach Science, New York, 1988). 
[10] T. B. Pittman, D. V. Strekalov, D. N. Klyshko, M. H. Rubin, A. V. Sergienko, and Y. H. Shih, Phys. Rev. A 53, 2804 
(1996). 

[11] M. D'Angelo and Y. H. Shih, quant-ph/0302146. 

[12] Recall that a set of objects S is said to be convex if for every pair of points a,b £ S and any t £ (0, 1) the points ta+ (1 —t)b 
are also in the set. That is to say, the straight line connecting any two point in the set lies in the set. 

[13] The completely random state has the property that probability of projecting on to an arbitrary state is independent of the 
state. In the case in the text, tr (pDo) = 1/N for any state p. 

[14] A. Sanpera, R. Tarrach and G. Vidal, Phys. Rev. A 58, 826 (1998), K. Zyczkowski, et. al, Phys. Rev. A 58, 883 (1998), 
A. Pittenger and M. H. Rubin, Phys. Rev. A 62, 032313 (2000). 

[15] K. Kraus, States, Effects, and Operations, ( Springer- Verlag, Berlin, 1983). 

[16] A. Peres, Quantum Theory: Concepts and Methods, (Kluwer Academic Publishers, Dordrecht, 1993). 

[17] N. Gisin, Helv. Phys. Acta 62, 363 (1989), L. P. Hughston, R. Jozsa, and W. K. Wooters, Phys. Lett. A 183, 14 (1993). 
[18] M. H. Rubin, Phys. Rev. A 61, 022311 (2000). 

[19] T. Aichele, A. I. Lvovsky, and S. Schiller, European Physical Journal D 18, 237 (2002). 

[20] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 

[21] B. M. Terhal, Phys. Lett. A 271, 319 (2000). 

[22] A. Pittenger and M. H. Rubin, Lin. Alg. Appl. 346 , 75 (2002). 

[23] O. Giihne, P. Hyllus, D. Bruss, A Ekert, M. Lewenstein,C. Macchiavello, A. Sanpera, Phys. Rev. A 66, 062305 (2002) and 

quant-ph/0205089. 
[24] A. O. Pittenger and M. H. Rubin, Phys. Rev. A 67, 012327 (2003). 
[25] M. H. Rubin, Phys. Rev. A 54, 5349 (1996). 

[26] J. W. Goodman, Introduction to Fourier Optics, (McGraw-Hill Publishing Co., New York, 1968). 

[27] A. Gatti, E. Brambilla, and L. A. Lugiato, Phys. Rev. Lett. 90, 133603 (2003) and quant-ph/0210109, 

[28] A. VanderLugt, Optical Signal Processing, (John Wiley & Sons, Inc., 1992). We use slightly different notation than that 

used in this text. 
[29] S. Popescu, Phys. Rev. Lett. 74, 2619 (1995). 



12 




FIG. 1. Illustration of an SPDC experiment. In arm A there is a lens LI a distance d± from the output face of the SPDC 
crystal and an aperture T a distance d[ behind it. A fixed detector Dl that collects all the light that is transmitted through T. 
In arm B the D2 is a planar array of detectors located a distance from the crystal. The coincidence circuit has been omited. 




FIG. 2. Illustration of a classical analog of the experiment in fig.l. The detectors are the same as those in fig. 1. The source 
is understood to be a stocastic source that emits pairs of identical pulses into A and B. The transverse wave number varies 
randomly from pulse pair to pulse pair. The lenses can be thought of as part of the source or as part of the optical systems 
in the two arms. However, it is important that T lies in the focal plane of the lens LI and the detector array plane lies in the 
focal plane of L2. 
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FIG. 3. Illustration of the ghost interference experiment. The aperture T is a double slit and the detector Dl is in the far 
field and only collects light with zero transverse wave number. 
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